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1 Introduction 



f(R) gravities as modified gravity theories [IH1] have much attentions as one of promising 
candidates for explaining the current and future accelerating phases in the evolution of 
universe [5]. It is known that f(R) gravities can be considered as general relativity (GR) 
with an additional scalar field. Explicitly, it was shown that the metric-/ (R) gravity is 
equivalent to the ojbd = Brans-Dicke theory with the potential, while the Palatini-/ (R) 
gravity is equivalent to the cubd = —3/2 Brans-Dicke theory with the potential [6]. Although 
the equivalence principle test (EPT) in the solar system imposes a strong constraint on f(R) 
gravities, they may not be automatically ruled out if the Chameleon mechanism is employed 
to work. It is shown that the EPT allows f(R) gravity models that are indistinguishable 
from the ACDM model (GR with positive cosmological constant) in the evolution of the 
universe [7]. However, this does not imply that there is no difference in the dynamics of 
perturbations [S|. 

On the other hand, the Schwarzschild-de Sitter black hole was obtained for a positively 
constant curvature scalar in [8] and other black hole solution was recently found for a non- 
constant curvature scalar j9j. A black hole solution was obtained from f(R) gravities by 
requiring the negative constant curvature scalar R = R [TO]. If 1 + f'(Ro) > 0, this black 
hole is similar to the Schwarzschild-AdS (SAdS) black hole. Also, its seems that there 
is no sizable difference in thermodynamic quantities between f(R) and SAdS black holes 
when using the Euclidean action approach and replacing the Newtonian constant G by 
G eS = G/(l + f'(R )). 

In order to obtain the constant curvature black hole solution from u f(R) gravity cou- 
pled to the matter", the trace of its stress-energy tensor T^ v should be zero. Hence, two 
candidates for the matter field are the Maxwell and Yang-Mills fields. Concerning the f(R)- 
Maxwell black hole, the authors [TU] have made an mistake to show the correct solution [TTj . 

In this work, we study the /(-R)-Maxwell black hole and its all thermodynamic quanti- 
ties, which are similar to the Reissner-Nordstrom-AdS (RNAdS) black hole when making 
appropriate replacements. We obtain the topological /(-R)-Maxwell black holes. Impor- 
tantly, we obtain the topological /(i?)-Yang-Mills black holes, which are similar to the 
topological Einstein- Yang-Mills (dyonic) black holes in AdS space. Since there is no ana- 
lytic black hole solution in the presence of Yang- Mills field, we obtain asymptotic solutions. 
Then, we confirm the presence of these solutions in a numerical way. 
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2 /(i?)-Maxwell black holes 



Let us first consider the action for f(R) gravity with Maxwell term in four dimensions 

1 



\fM 



16irG 



tfx^g R + f(R) - F^F^ . (2.1 



From the variation of the above action (12.11) . the Einstein equation of motion for the metric 
can be written by 

R^(l + f'(R)) -\(R+ f{R))gp, + (<7^V 2 - V,vJ)f'(R) = 2T^ (2.2) 
with the stress-energy tensor 

T^ = F w F v ?- g -^F pc FP° withT% = 0. (2.3) 
On the other hand, the Maxwell equation takes the form 

V M F^ = 0. (2.4) 
Considering the constant curvature scalar R = Rq, the trace of (12. 2 p leads to 

Ro(l + f'(Ro)) -2(r + f(Ro)^ = (2.5) 
which determines the negative constant curvature scalar as 

*-7lpT s4A ' <a (26) 

Substituting this expression into (12. 2p leads to the Ricci tensor 

2 

which implies that R^ v ^ ^fd^ (pure AdS 4 space) unless T^ v = 0. 
We introduce a static spherically symmetric metric ansatz, 

dv 

ds 2 = -N(r)dt 2 + + r 2 dn 2 2 (2.8) 



and a gauge field as a solution to ( 12. 4p 



Mr) = — ~ - (2-9) 

r + r 



which provides an electrically charged black hole with A t (r + ) = 0. Solving the Einstein 
equation ( 12. 2 p together with the condition of constant curvature scalar, we obtain the 
solution for a metric function 

2GM Q 2 R 9 , 

N( r ) = i 1 -5 V. (2.10) 

U r + (1 + !' (R ))r 2 12 1 UJ 

We note that the topological /(-R)-Maxwell black hole solution is also found to be 



2GM t Q 2 _R 

1 + f'(R ))r 2 ~ VI 



N k (r) = k + - ■ ^— . - - -V (2.11) 



when considering the metric ansatz 

ds 2 = N k (r)dt 2 + N;\r)dr 2 + r 2 dZ 2 , (2.12) 

with dT? k = d6 2 + a^.(9)dip 2 . Here (7^(0) denotes sin^, 9 and sinh^ for k = 1 (spherical 
horizon), (flat horizon), and k = —1 (hyperbolic horizon), respectively. 

We could derive all thermodynamic quantities since the analytic solution was known as 
(12.1 Op . First of all, the Hawking temperature is calculated to be 



N' 1 
T H (r + ,Q) = —\ r ^ r+ = — 



Q 2 Ror+ 



(2.13) 



In order to compute other thermodynamic quantities, it would be better to use the Eu- 
clidean action approach [13] because we are working with f(R) gravities. To make the 
action Euclidean, the time coordinate should be made imaginary by substituting t = it. 
In this case, to eliminate the conical singularity at the horizon r = r + , the coordinate r 
should be periodic with the period j3 = l/T H . For this purpose, we have to calculate the 
Euclidean action [I~5l[T6] 

AS? = Sf M + S GH + S ct + S cFl (2.14) 
where the Euclidean bulk action is 

Sf M = ~i^~Q J #Xy/B[R + f(R) - F»»F" U ] ■ (2-15) 

Here Sqh is the Gibbons-Hawking term to make the variation at the boundary clear and 
S ct is the counter term for asymptotic AdS4 space. We are working with the canonical 
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ensemble as the fixed charge ensemble. In this case, we need to introduce the charge-fixing 
(cF) as a boundary surface term [TS] 

1 



S, 



cF 



AtxG 



d 3 xVhF» v n v A v 



(2.16) 



where hij is the induced metric on the boundary surface and is a radial unit vector 
pointing outwards. If one does not introduce S c f, one is working with the grand canonical 
ensemble. Also, the extremal black hole whose horizon is degenerate is considered to be the 
ground state in the canonical ensemble [T7]. The location r + = r e of extremal horizon is 
determined by the condition of Tjy(r e , Q) = 0. That is, in order to derive the Helmholtz free 
energy, we must subtract the extremal mass M f = Ml from fl2~T4]) . Taking into account 
all leads to 



ASf - Ml 



0{l + f{Ro)) 



48G 



12r 



36Q 2 



;i + f'(Ro))r, 



24 r 



5 
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(2.17) 
(2.18) 



(l + f'(R ))r e 
= = 0Ef - S f BH . 

Here is the Helmhotz free energy, (3 is the inverse of the Hawking temperature, and S^ BH 
is the Bekenstein-Hawking entropy. The energy, Bekenstein-Hawking entropy, and heat 
capacity are given as 

8(AS E - Ml) 



E f (r + ,Q) 



df3 

;i + fgoK 

2G 



M f {r + ,Q)-Ml 

Q 2 



of 

J R 



BH 



C f (r + ,Q) 



PEf - PFl = (1 + f(R )) 

fdEf\ 
V&TJq 

2(1 + rWjTTTj. 

G 



rmy + 

A(r+) 



2 

12 + 



-Ml, 



AG 



4r 'i + TT7TO + R °< 



(2.19) 
(2.20) 



4ri - 



12Q 2 



+ Ror 4 _ 



where 



M/ = M/(r e ,Q) = i±^ 



2Q 2 



r e (l + f'(R )) 



(2.21) 



(2.22) 
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Figure 1: Thermodynamic quantities of the /(-R)-Maxwell black hole as function of horizon 
radius r + with fixed Qf = 1, £ = 10, and G e ff = 1 ; temperature T H , heat capacity C, and 
Helmhotz free energy F . 

is the mass of the extremal black hole and A(r + ) = Atir^ is the horizon area. In this case, 
E measures the energy above the ground state. 
Considering replacements of 



G 



12 

G eff , R = 4A f -— 



(2.23) 



l + f'(Ro) u J £ 2 ' l + /'(i2o) 

the f(R)- Maxwell black hole becomes the RNAdS black hole exactly. In this case, the 
ADM mass , Hawking temperature Th, and the Bekenstein-Hawking entropy S BH take 
compact forms 



M f (r + ,Q f ) 



2G, 



off 



1 + — + — 

rl ^ P 



T H (r + ,Q f ) = - 



Q 



3r + 



of 
°BH 



7rn 



Finally, the heat capacity C* and Helmholtz free energy 

F f 

are given by 



C f (r + ,Q f ) 
F f (r + ,Q f ) 



2-nr\ 
1 



3r 4 + + £ 2 (r 2 + -Q 



3r\ + £ 2 (-r< 



3Q 2 



4G cff r + 



r 2 , +3Q 



t 2 



f 
- M P 



G c s 
(2.24) 



(2.25) 



(2.26) 



At this stage, we have to mention the other thermodynamic quantities obtained directly 
from the metric function N(r) in (12.101) . In this case, all thermodynamic quantities of 
M, Sbh, C, F are obtained from the replacements as 



12 

°0 — > mi 



1 + f'(Ro) ^ Q f 



(2.27) 



There exists a slight difference in M, Sbh, C, F between G in the direct method and G e fj 
in the Euclidean action approach. The first law of thermodynamics is satisfied for both 
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cases 



dM f = T H S f BH , dM = T H S BH . (2.28 
One curious quantity derived from f(R) gravities is the Bekenstein-Hawking entropy 

of 



BH 



) (2-29) 



which was also derived from the Wald method [6]. On the other hand, the conventional 
Bekenstein-Hawking entropy is 



nr 2 



Sbh = (2.30) 

For example, if one uses S BH to check the first law of thermodynamics, one immediately 
finds that it is not satisfied as follows 

dM ^ T H S f BH . (2.31) 

At this stage, there is no way to test which approach provides the correct thermodynamic 
quantities for /(i?)- Maxwell black holes. Anyway, the entropy issue should be resolved. 

The global features of thermodynamic quantities are shown in Fig. 1 for G c s = 1 = G. 
Under this setting, there is no difference between two approaches: M$ = M, S BH = 
Sbh, = C, = F. From the first and second graphs, we observe the local minimum 
T H = T (C blows up) at r + = r , in addition to the zero temperature T H — (C = 0) at 
the extremal point of r + = r e and the maximum value Tjj = T m (C blows up) at r + = r m 
known as the Davies point. We note a sequence of r e < r m < vq. For r e < r + < r m , the 
black hole is locally stable because of C > 0, while for r m < r + < r it is locally unstable 
(C < 0). For r + > r , the black hole becomes stable because of C > 0. Based on the local 
stability, the /(i?)- Maxwell black holes are split into small black hole (SBH) with C > 
being in the region of r e < r + < r m , intermediate black hole (IBH) with C < in the region 
of r m < r + < ro, and large black hole (LBH) with C > in the region of r + > vq. 

Importantly, the free energy from the last graph in Fig. 1 plays a crucial role to test 
the phase transition. A black hole is globally stable when C > and F < 0. We note that 
F = at r + = r e , because of F = M — M e — ThSbh with Tn(r e , Qf) — 0. We observe two 
extremal points for free energy: the local minimum F = F m i n at r + = r m and the maximum 
value F = F max at r + = r . The free energy is negative for r e < r + < r m and it increases in 
the region of r m < r + < r . For a point of r + = r\ > r , it is zero and remains negative for 
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r + > r\. The temperature of T — T\ (determined from the condition of F = 0) at r + = ri 
may play a role of the critical temperature in Hawking-Page phase transitions. The related 
phase transition was discussed in Ref. [18]. It can be shown that the Hawking-Page phase 
transition II between SBH and LBH unlikely occurs in the /(i?)-Maxwell black holes. 

3 /(i?)-Yang-Mills black holes 

We consider the action of f{R) gravity coupled to SU(2) Yang-Mills field in four dimensions 

SfYM = ^ J d'x^g {R + f(R) - F; V F^} , (3.1) 

where F" = d^A® — d v A a ^ + e abc A b A c v . From the action (jXTJ , the Einstein equation of 
motion can be written by 

iV(l + f(R)) ~ \(R + f(R))g,u + (g»uV 2 - V,V,)f (i?) = 2T™ (3.2) 

with T^, M the stress-energy tensor for the Yang-Mills field. For the constant curvature 
scalar R = Rq, taking the trace of (I3.2p leads to 

Ro(l + f(Ro^ - 2(i? + /(i2o)) = (3.3) 

which determines the constant curvature scalar as 



2f(Ro) 

fW - 1 



Rq = £ y^> . = 4A/ < o. (3.4) 



Now we consider the topological metric ansatz 

dsl = -e 2<t>{r) N(r)dt 2 + N~\r)dr 2 + r 2 dE 2 k , (3.5) 

with dT^l = d9 2 + o\{ff)dif> 2 . A dyonic solution ansatz for Yang-Mills gauge field is given by 

A = {u{r)r 3 dt + u(r)nde + d e a(9)r 3 + a(9)cu(r)r 2 dcp}, (3.6) 

where u(r)[u(r)] describe the electric [magnetic] charged configurations and is the Pauli 
spin matrices for SU(2). 
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Substituting ( 13. 5 j) and ( 13. 6 p into the action (13. ip . and after variations with respect to 
N, 0, cu, u, one finds their equations of motion: 



S N S; 2r(l + /'(iZ))^ - r 2 {/'"( J R)( J R'(r)) 2 + f"(R)(-cf>'R'(r) + iT(r)) 



- 4 (-') 2 j^— = 0, 

2(w 2 - kf 



5^S; -2k + 2rN' + 2N- r 2 f(R) + AN(oo') 2 + {U ^ ] + 2e" 2 * (r V) 2 



(3.7) 

2w 2 c^ 



iV 



+2r 2 Nf"(R)(R'(r)f + r/^^^M'lr) + 2rNR"(r) + rN'R'(r) 
+rf'(R) \ - 2rN((f)') 2 - ANcj)' - rN" - 2rN(f)" - 2N' - ZrN'<j>'\ = 0, 



6 W S; r 2 Noo" + r 2 (N' + 0'iVy - w(w 2 - k) + 



e 2cl, r 2 u 2 uj 



5 U ,S; rV + (-r 2 0' + 2r) «' 



2w 2 u 
TV 



(3.8) 
(3.9) 
(3.10) 



where the curvature scalar R(r) is given by 



R{r) 



2 + rN'(A + 3r0') + r 2 iV" + N{2 + 2r 2 (0') + 4r0' + 2r 2 (j>") 



(3-11) 



Note that the prime (') in f(R) and N,u,<j) denotes the differentiation with respect to 
R and r, respectively. It is a formidable task to solve the above four equations directly. 
Therefore, we consider the constant curvature scalar which implies that 



R(r) = R 0l R'(r) = R"(r) = 0. 



(3.12) 



Actually, we have used the condition (I3.12p to derive the /(i?)-Maxwell black holes in the 
previous section. Plugging (I3.12p into the four equations leads to simplified equations 



5 N S; r(l + f( J R o ))0 / -2i 



J\2 



2e- 2 ^u 2 uj 2 



N 2 



0. 



S+S; -2k + 2rN' + 2N -r 2 f(R ) +AN(u 



2{u 2 -k) 2 _ _ nJ , ( , , NO 2u l u 



+ 2e 



(3.13) 

2, ,2' 



20 / r 2 {v!) 2 + 



N 



+rf'(R ) ( - 2rN((j)') 2 - ANcj)' - rN" - 2rN(j)" - 2N' - 3rN'(j)') = 0, (3.14) 
6 U S; r 2 Nu" + r 2 (N' + <f>'N)u,' - cu(cu 2 - k) + = 0, (3.15) 



N 



6 U S; r 2 u" + (-r 2 0' + 2r) u -2 



<jJ 2 U 

~N~ 



(3.16) 



Note that for f(Ro) = /'(-Ro) — 0, they reduce to those derived from the topological 
Einstein- Yang-Mills theory [2T]. It is well known that there is no analytic black hole solution 
to the Einstein- Yang-Mills theory. Hence, we can find either the asymptotic solution with 
finite terms or the numerical solution. 

First, we wish to derive the asymptotic solution at infinity of r — > oo. Equations (13. 13ft - 
( I3.16P can be solved by considering asymptotic forms for metric and gauge field functions 
up to A-order 



k 



N(r 
m(r 



2m(r) 



Ro 2 
12 ' 



Mi M 2 M 3 
M + — + + ^ 



r 
r 



M i + m + 



0J 2 



0J 3 OJa 



UJ 5 



u{r) = u 00 ^ ^^" + -7 + ^- + ^ + ° 

iy ry £ ty<J ly^i ryi) 



(3.17) 
(3.18) 

(3.19) 

(3.20) 



where M, Uoo, u\, Uoo, and u\ are five constants evaluated at infinity and other Mi, Ui, 
and Ui are expressed in terms of these constants and 1 + /'(-Ro) appeared in Appendix A. 
Let us compare /(_R)-Yang-Mills (fYM) black holes with Einstein- Yang-Mills (EYM) black 
holes. We observe the relations of coefficients between two black holes 

M EYM 



M, 



fYM 



fYM 



i + f'(R 

EYM 



UJ 



u7 M 



for i = 1,2,3,4 
M ? YM for^ 



2,3,4. 



(3.21) 
(3.22) 



It seems that there is no longer simple relations between two black holes for i > 5. Hence, 
it is not easy to derive any concrete form for thermodynamic quantities of /(i?)-Yang-Mills 
black holes. Exceptionally, the form of Hawking temperature can be derived to be 



T; 



fYM 
H 



k 



4 



2 



2m' (r j 



(3.23) 



because it will be determined by the variables defined at horizon. Using (14. lip , it takes the 
form 

1 



T, 



fYM 
H 



, Ro 2 (k-ul) 2 
k -rl 



4- ri(l + /'(i2o)) • (3 ' 24) 

In the case of purely magnetic charged black hole with ii = and /'(-Ro) = 0, it reduces to 
Eq.(28) in Ref. [21]. Here uq = u'{r + ) may be considered as a counterpart of A' t {r + ) = Qjr\ 
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in the /(-R)-Maxwell black holes. Furthermore, one finds the metric function N (m(r) 
M + Mi/r) up to -^-order and gauge field functions u and u up to —order 



2M 


JQ 2 M + 


r 




J 




+ -, 


u(r) ~ u, 


r 





cj ~ Woo + — , w(r) ~ Uoo - — (3.26) 

r 



with the Yang-Mills magnetic charge Qm = k—u^ [21] and the Yang-Mills electric charge Q. 
Here we reset U\ = J and U\ = —Q to make a connection to holographic super-conducting 
models using the AdS/CFT correspondence [22] and [23] for higher dimensional cases. 
Using the holographic interpretation with = 0, Uoo is the chemical potential, Q is the 
electric charge, and J is the component of the current Jj on the boundary at infinity which 
is connected with the spontaneously broken part of the bulk gauge symmetry. 

We note that ( I3.25P and (I3.26P with Q 2 M = and = Q/r + [imposed by u(r + ) = 0] 
reduce to those of the topological /(i?)-Maxwell black holes when turning off the mag- 
netic charge gauge potential and setting Ui = 0(i > 2). Especially, the constant u\ = J 
corresponds to an order parameter describing the deviation from the Abelian solution of 
/(#)-Maxwell black holes. 

Finally, it is also interesting to explore the other case of purely magnetic charged black 
holes obtained by choosing k = 1, u(r) = 0. Its asymptotic solution appeared in Appendix 
B. In the case of Einstein- Yang-Mills black holes, these black holes are stable against grav- 
itational and sphaleronic perturbations for tu + > l/y/3 = 0.577 for large |A| [TH]. Actually, 
the stability condition corresponds to that a gauge field u(r) has no zero. Hence, we conjec- 
ture that purely magnetic charged black holes in the f(R)- Yang- Mills theory has a similar 
property because for constant curvature scalar and 1 + f'(Ro) > (no ghost condition), the 
/(-R)-modification to the Einstein- Yang- Mills black hole will be minimized. We will check 
in the next section that the zero of oj(r) appears only for u + < 1/ v3 = 0.577. Furthermore, 
there exist nodeless solutions for k — 0, 1 Einstein- Yang-Mills black holes (21], which means 
that these black holes are stable. 

In the next section, we will find numerical solutions for k = 1, dyonic black holes and 
k — 1 purely magnetic charged black holes. 
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4 Numerical results 



We numerically solve fl3.13p - fl3.16p with boundary conditions of (I3.25P and (I3.26P using a 
standard shooting method in Mathematica®7 . The k = 1 Einstein- Yang-Mills black holes 
was discussed in Ref. [ISIED], while k = 0,-1 Einstein- Yang-Mills black holes was found 
numerically in Ref. |21j . 

In order to obtain numerical solutions, we transform equations fl3.13p -f l3.16p into 

2L 2 u 2 e- 2 * + (u') 2 N 2 ) 



m'(r) 



2r 2 (u 2 u 2 e- 2 ^ + {u') 2 N 2 ^j +N^{k- u 2 ) 2 + r 4 ((u e -*)') 
2r 2 (l + f'(Ro))N 
2rW (cu 2 u 2 e- 2 't > + (cu') 2 N 2 ^) 2u 2 (We" 2 ^ + (u/) 2 iV 2 ) ' 



;i + f(i? )) 2 iv 2 (i + /' (i? )) 3 iv 4 



(4.2) 



,/ / N'u' cu(k-cu 2 ) uu 2 e- 2 * , A , 

„ >H0V -^ + ^, (4.4) 

where we used the relation H3.4[) to include 1 + f'(Ro) only as /(i?)-gravity effects. First, 
let us develop the solution forms near the non-degenerate horizon at r = r + . Because of 
N(r + ) = 0, we derive a relation of u(r + )u(r + ) = from Eq. (I4.4p . Choosing u(r + ) = 0, 
it is easily shown that w'{r + ) = u"(r + ) = 0, which implies that u(r) = 0. This is not the 
case. So we choose u{r + ) = instead, and the solution near the non-degenerate horizon 
can be expanded as 

0(r) = </>+ + 0'(r+)(r - r+) + 0(r - r+) 2 , (4.5) 

m{r) = m + + m'(r + )(r — r + ) + 0{r — r + ) 2 , (4-6) 

u{r) =u+ + u\r + )(r - r+) + 0(r - r + ) 2 , (4.7) 
u"(r ) 

u(r) = u (r - r+) H \r^-{r ~ r +) 2 + °( r ~ r +) 3 ' ( 4 - 8 ) 
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where the coefficients are determined by equations 



m+ = m(r+) = — ( 1 - —r + ) , (4.9) 
2ul Uk - u\y + r%ute 

0/(r+H MTWT - (410) 

m'fr ) - ( fc ~ ^+) 2 + <^' 20+ f4 n) 

«'(r + ) = "-^f , (4-12) 
tt »f r ) - _^ f i _ rW + (4k - R rl) \ 

which satisfy at the horizon r = r + . Since the metric function is zero at the horizon 
r = r + and it should be positive outside the horizon, we have to choose a condition of 
N'(r + ) > 0. This restricts the range of u + through the inequality 

2m'(r+) < k - (4.14) 

which yields a positiveness of u/(r+) > for u; + ((X>^ — k) > 0, while m'(r + ) > for 
1 + /'(-Ro) > 0. Note that 0'(r + ), m'(r + ), u/(r+), and w"(r + ) depend on r+, 0+, o;+, and 
Mo, which means that they are four independent parameters describing the near horizon 
geometry of (I4.5p — (14. 8p . On the other hand, there are five independent parameters of 
M, Uoo, u>i, Uoq, and u\ describing asymptotic region of ( I3.18p -(l 3.20p . Remembering that 
(I4.ip -( |4T1| are two first- and two second-order differential equations, we need six initial 
parameters to solve the equations numerically at each boundary of horizon and asymptotic 
infinity. However, considering N(r + ) = 0, we choose u(r + ) = and then, (14. 3p leads to 
a first-order differential equation. This is why four independent parameters is enough to 
specify the near horizon geometry of /(i?)-Yang-Mills black hole. In addition, we have 
time-rescaling symmetry so that we can replace <fi by + O without loss of generality. This 
means that either + or 0(oo) can be cast to zero. Here we choose 0(oo) = to achieve 
an asymptotic AdS4 space which is similar to that of /(i?)- Maxwell black hole solution. 
Hence, the asymptotic solution has five parameters less than six of general analysis. In this 
manner, we show that <p + is not an arbitrary parameter. 

We are now in a position to solve the initial value problem, for given r + , _R , uj + , and 
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(a) k = +1, uj + = 1.08, and u = 0.33 (b) k = 0, lu + = 0.1, and m = 0.77 

Figure 2: The numerically solved functions m (solid lines), u (dashed lines), u (dotted 
lines), and (dot-dashed lines) are depicted with respect to radius log 10 r for r+ = 1 and 



R = -3.6 (a = -0.9, = 1.125, c x = 3.025). [(a)] For fc = 1, we have m + = 0.65 and 
+ = —0.23 with 0(oo) = 0. Then, we obtain = 1.0 so that Qm = (1 — = 0. (b) 
For fc = 0, we get m + = 0.15 and 0+ = —0.11 with 0(oo) = 0. Then, we obtain Uoo = so 
that Qm = 0. 

Uq, by introducing a specific form of f(R) gravity as 



R 



f(R) = (4.15) 



proposed in Ref. [24] and setting n — 1 for simplicity in this work. Imposing the constant 
curvature scalar (]3.4p . Eq. fl4.15p implies 

Cl = \+md ■ (416) 

The Einstein limit exists in the constant curvature case, which shows that c\ — > 2/3 Af /a as 
a — > 0. Then, f(R) becomes a negative cosmological constant (f(R) — > —2Af). In what 
follows, however, we consider only a nonvanishing a case. 

The numerical solutions to ( I4.ip - (l4.4p are depicted in Fig. |2] graphically for a = —0.9, 
(3 = 1.125, Ci = 3.025, for which Rq and f'(Ro) are fixed to either Rq = —3.6 and f'(Ro) = 
—0.1 or R = 1.8 and f'(Ro) = —10. In this work, since we are interested in asymptotically 
AdS space with 1 + f'(Ro) > 0, we choose the former of Rq = —3.6 and f'(Ro) = —0.1. 
Choosing the horizon radius to be r + = 1, the left and right figures are distinguished by 
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0.5 - 



1 2 3 

(a) k = +1, uj + = 1.08, and + = 



— l°gu, r 
-0.03 



1.2 














, 




m(r) 


1.0 










0.8 


/ 






e m 


0.6 










0.4 










0.2 


w(r) 













1 


2 


3 



lo g| r 



(b) k = +1, u> + = 0.1, and <f) + 



-0.01 



Figure 3: The numerically solved functions m (solid lines), u (dashed lines), and (dot- 
dashed lines) of purely magnetic case (u(r) = 0) are depicted with respect to radius log 10 r 
for r + = 1, and Ro = —3.6 (a = —0.9, (3 = 1.125, c\ = 3.025). For convenience, k — 1 case 
is plotted only, and we have m + = 0.65. 



specifying remaining parameters k, u + and uq\ (a) k = 1, lo+ = 1.08 and uq = 0.33 (b) 
k = 0, u + = 0.1 and uq = 0.77. Furthermore, these numerical solutions are being used to 
determine the form of parameters in asymptotic AdS space. Matching its asymptotic forms 
( 13723]) and flQgp , we find that for [(aJJ M »s 4.07, w 1.00, wi w 2.63, w 1.38, and 
Wi w -4.14, while forgbJJ M w 0.58, w 0.00, w x w 0.13, w 0.87, and m w -0.87. 
At this stage, we point out that that the magnetic charge Qm = {k — oo^) vanishes for both 
cases, so that their asymptotic geometry are similar to the /(i?)-Maxwell black holes. 

It is also interesting to explore the other solution of purely magnetic charged black holes 
numerically by choosing k = 1, u(r) = 0. Considering its asymptotic solution appeared in 
Appendix B, we find the numerical solutions. For a given /(i?)-form (I4.15p . the numerical 
solutions to Eqs.(5.2)-(5.4) could be developed for the same values as in the dyonic black 
hole solution. Setting the horizon at r + = 1, the two graphs in Fig. 3 are distinguished 
by specifying a remaining parameter uj + : 



1.08 (b) 



0J A 



0.1. Furthermore, this 



numerical solutions are used to find the parameters in the asymptotic solutions. We find 

-2.54 for 



M w 1.19, Woo ~ 



1.17, ui 



while M w 1.21, co^ w -0.02, u x w 0.14 for 
(b) The stability condition for this black hole corresponds to the condition that a gauge 
field u(r) has no zero. As is shown in the Fig. 3, we find that the zero of u(r) appears 
for tu + < l/VS = 0.577. Hence, we conjecture that the stability condition for the Einstein- 
Yang- Mills black hole holds for the /(i?)-Yang-Mills black holes with 1 + f'(R ) > 0. 
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Table 1: Asymptotic forms of numerical solution for f(R) -Yang- Mills black holes 



dyonic solution 



k = +1 
(u+ = 1.08, n = 0.33) 



k = 
(u+ = 0.1, n = 0.77) 



purely magnetic solution 



oo + = 1.08 
(* = +l) 



uj + = 0.1 
(k = +l) 



m(r) 
uj{r) 
u(r) 



4.07- 



1.00 



18.87 

r 

2.63 



1.38- 



r 

4.14 



r 

15.57 

r 4 



0.58- 



0.00 



0.87- 



0.42 

r 

0.13 

r 

0.87 



0.01 



1.19 



1.71 



4.18 

r 
2.54 

r 





3.58 



1.21 - 



-0.02 



0.56 

r 

0.14 

r 





0.01 



5 Discussions 

First of all, we summarize all numerical solutions to /(i?)-Yang-Mill black holes in Table 1. 
This table shows asymptotic solution forms constructed in the numerical way: two dyonic 
solutions for k = 1 and k = black holes and two magnetically charged black holes for 
k = 1 and u + = 1.08, 0.1. The former was developed to compare with the topological 
Einstein-Maxwell black holes, while the latter was displayed to see the stability of f(R)- 
Yang-Mills black holes. We find that for 1 + /'(_R ) > 0, the /(i?)-Yang-Mills black holes are 
similar to Einstein- Yang-Mills black holes in AdS space. In this case, one may develop the 
second-order phase transition between /(i?)-Maxwell and f(R)- Yang- Mills black holes to 
explain the holographic superconductor without Higgs field as in the Einstein theory |22j. 
The difference is that the AdS4 space was constructed not by introducing a cosmological 
constant, but by choosing an appropriate f(R) function in (13.41) . Also, it seems that for 
1 + f'(Ro) > 0, the stability condition of magnetically charged Einstein- Yang-Mills black 
holes holds for magnetically charged f{R)- Yang-Mills black holes. 

The condition of 1 + f'(Ro) > is related to no ghost state for graviton propagations 
on AdS4 space [251126] . the positiveness of effective Newton constant G c ff > in cosmo- 
logical implications [6],[7], and a necessary condition that f(R) black hole becomes a type 
of Schwarzschild-AdS black hole [TU]. In this work, this condition is necessary to obtain 
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/(-R)-Maxwell black hole and to derive its thermodynamic quantities. Also, /(i?)- Yang- 
Mills black holes requires this condition to have asymptotic and numerical solutions. The 
other condition of f"(Ro) < is not necessary to obtain the constant curvature black hole 
solutions. However, this condition may be needed to be free from the Dolgov-Kawasaki 
instability related to tachyonic mass [21 [6] in the perturbation analysis of /(i?)-Maxwell 
(Yang-Mills) black holes. We hope to make a progress on the perturbation analysis. 

At this stage, we would like to mention the close connection between f(R) and Einstein 
black holes by rewriting the action (12. ip as 



{R+f(R)}-F^F 



16nG 

In this case, Einstein equation takes the form instead of (I2.7P 



(5.1) 



R/AV ~ A /fiV + l + flfJl\TlW ( 5 ' 2 ) 

with T^y = 47^. Introducing a replacement of G e ff —> Gj (1 + f'(R )), the above equation 
become 

Rpv = A f g^ + STrGefff^. (5.3) 

The solution is determined by 



H 12 



which is the same form as (I2.10p in the unit lQnG = 1. Also, we may make such a 
replacement for /(i?)-Yang Mills theory by rewriting ( 13. ip as Sjym- In this case, the 
Mi, Ui, and it, in Appendix A including the substitution rules ( I3.2ip and ( I3.22p may be 
conjectured by following v = 4-nG/e 2 with e 2 = 1 + f'(Ro) in Ref. poj, where — F 2 /4 was 
used instead of —F 2 . 

Finally, we wish to comment on two points. One is to answer to the question of "is it 
possible to apply the reconstruction technique of Ref. [27] which was developed for pure 
f(R) gravity to the f(R) with Maxwell (Yang-Mills) field?". The answer is "yes" because 
the pure metric f(R) gravity is equivalent to the o>bd Brans-Dicke theory with the potential 
term. Expressing = 1 + f'(R), we transform (12. ip and ( 13. ip into the Brans-Dicke theory 
[0, V{4>) = f — Rf'(R)] with Maxwell (Yang-Mills) field. As far as the constant curvature 
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scalar black hole solution is concerned, we can obtain the same black hole solution from 
the reconstructed Brans-Dicke theory with Maxwell (Yang-Mills) field [251121] • The other 
is to answer to the question of "in the case of coupled YM- f(R) theory |30j, 



b~g l In 



—0.5F a F^ 0, 



Symj = j d 4 *V=9 |^ - 4ttG(1 + f{R))F; v F^ 

■(5.5) 

do we expect to obtain the similar solution?". In the case of f(R) = 0, the last term of 
the above action reduces to the effective Lagrangian of SU(N) Yang-Mills theory up to 
one-loop order with 

" = k&T N - (5 ' 6) 

Even though its equation of motions take complicated forms, we expect to have similar 
numerical solution found here for the constant curvature scalar black hole. This may 
be true because for the constant curvature scalar black hole, the replacement of G e ff — > 
G/(l + f'(Ro)) is expected to make the numerical solution simple unless b log-term plays 
an important role. 

Consequently, the /(-R)-MaxweH (Yang-Mills) black holes imposed by constant curva- 
ture scalar and 1 + f'(Ro) > are closely related to the Einstein-Maxwell (Yang-Mills) 
black holes in AdS space. 
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Appendix A: Coefficients for asymptotic solution to 
/(i?)-Yang-Mills black holes 



Mi 
M 2 
M 3 
M 4 



k 2 + u\ - oo 2 R /G - 2{k + 12u 2 J Rq)^ + w, 



1 

CO 



2(1 + /'( J R )) 



1 + f (#o) 

le^MooWiW^ + i2oo;?(fc - 2co 2 J + ^{u 2 + (k - co 2 J 2 - 6*4 (fc + 2co 2 J/R ) 



Ro{l + f'{Ro)) 



1 



2(1 + f'(Ro)) [ 



M 



iU\U c 
Rn 



2 48M 00 ^{3Mu 00 -4ui(fc-6«y j Ro-2a;^)} 

Wl + ^ 



+2u^u; c 



16o;ia; 00 (2A ; 2 + ^-5^ + 3^) , 96^00^ (2k + u&) 



1 



30(1 + f'(Ro)) 



-Gut 



Ro 

12u) 2 (21k 2 + Gu\ - 99A&& + 78w 



Rl 
4 ) 

ooJ 



Rn 



576w 1 2 m^ (A; - 2u4) 288wiu; 0O wiw 0O (9£; - 48w 2 J J R - 28w 



+ 



Rl 



360w 1 w oo M(A; - A320u 2 o Mu 1 u oo + UAu^Uk - 12a&)(fc 



2 \2 



R 







Rl 



144a& M j(13fc-12a&) _ 17280^^ 41472^^ (2/c + 3a&) 



Rl 



Rl 



1728u 2 00 u 2 00 (k-u 2 J(13k + 6u 



2 ) 



it* 



uj: 



12{u\ + k 2 - 2fca4 + u, 
Ro 



1 ) 

CO I 



60(1 + /' (R )Y 
Rl 



J. + ^oo , 

Kq 
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6 

—ujoo(k - UuJRq - u4) 

Ho 

6 

— (-SmiMooWoo/^q + - 4u^/i2o - w£>)) 
K 



2R 

+6co 2 co 
1 



4(3M + UuiUoo/Rofai 
288u OD u 2 O0 {5k - 4w, 



12w 00 (7fc 2 - 2m? + l44M^/i2o - 10fca& + 3w, 



1 ) 



Rn 



2 ) 

DO/ 



6^ 



R>2 

576w 0O w 1 w 0O (llA; - 24u^/i? - 7u4) 144 Woo M(4A; - 60*4/i2o - 4w 



Rl 



Rn 



12wi {39fc 2 - + 144t4/i$ - + 27o4 - 264?4(fc - 2u 2 oc )/R } 



10(1 + f'(R )) 



24u\ 

Rn 



3ut 



Ro 

24k 2 - (48k + 432u 2 jR )u 2 oo + 24a;, 

Rn 



"2 
"3 
U4 

U5 



12 2 
/t 



ito 
2 

Rn 



ZuiUiUoo — 2 h u 00 u 1 



24u 00 a&(-fc + 



Rn 



5Ro 
4 



[48^^ {3M Woo + wi(-6fc + 24n^/i? + 7u&)}] 



30(1 + /'(i2o)) 
144a& {-3t& + (-6* + 60^/^ + 4a&)(l + f'(R ))Y 



R 2 
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Appendix B: Asymptotic solution to a magnetically 
charged black hole for /(/^-Yang-Mills theory 

In this case, SU(2) Yang-Mills gauge field is given by 

A = {u{r)ndB + [d g a(9)T 3 + a(9)cu(r)T 2 } d V } . (5.7) 

Three equations of motion for N, <p and u are given by 

6 N S- r{\ + f\R )W - 2(u/) 2 = 0, (5.8) 
8+S\ -2 + 2rN' + 2N- r 2 f{R ) + 4N{u') 2 + ^ ~ ^ 

+rf'(R ) ( - 2rN((f)') 2 - ANcj)' - rN" - 2rN(f)" - 2N' - 3rN'<j)'^ = $.9) 

5 U S; r 2 Nuj" + r 2 (N' + (j) 1 N)J - u(u 2 - 1) = 0. (5.10) 

The above equations can be solved by assuming asymptotic forms up to ^--order 

2m R n . „ „ . 

N = l- — -^r 2 , (5.11) 

Mi M 2 M 3 M 4 M 5 / 1 \ 

m = M H 1 — H — H H r + CM — , (5.12) 

^•J A i^jO i^t: ^0 V y^lO / 

uj^uj^^ 1 — 5-H — s-H — — =- + OI -B I , (5.13) 

ty ty£i tyO iy>^ iy<J \ /y\) J 
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where Wqo, oji, M are three constants and the coefficients Mj and Ui are expressed in terms 
of Woo, ui, M, and 1 + f'(Ro) 

M = ^ ~ 1)2 ~ ^° 
2(l + f(i? )) ' 



<2 



i + f'(Ro) 



J2o(l + /'(iZo)) 
-w?(M + 2u 1 u oo )R + 16^1^(1 - u£)(2 - 3a&) 



2i2o(l + /' (i? )) 



M 5 = - 5fi 2 (1 + fW ) W^i ~ SOM^i^C-l + a4) + 24(-13 + 12^)^00 - u^f 

-6R u 2 (7 - 33a& + 2&4)] + — - + ^ /(jRq))2 {-^i + 12(-1 + a&) a } 
6^oo(l - £«&) 

0^2 — — 



2 ) 



UJ4 



Rq 

6ui(1 — oj, 
Rq 

-3(2M + oj^^Rq + 6^(1 - u&)(7 - 3w 2 



CJ 3 

'to 



/? 2 



w 5 = [RoVi-WMu) 00 (-l+«% ) -6^(13-22^ + 9^ 



5i? 2 



{i? c 2 -8(-l + ^) 2 }. (5.14) 



10^0(1 + f'(Ro)) 
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